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ABSTRACT 

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena 
over timescales of seconds to (at least) hours, including nulling, subpulse drifting, 
and mode-changing. Studying pulsars like PSR B0823+26 provides further insight 
into the relationship between these various emission phenomena and what they might 
teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that 
PSR B0823+26 has a weak and sporadically emitting ‘quiet’ (Q) emission mode that 
is over 100 times weaker (on average) and has a nulling fraction forty-times greater 
than that of the more regularly-emitting ‘bright’ (B) mode. Previously, the pulsar 
has been undetected in the Q-mode, and was assumed to be nulling continuously. 
PSR B0823+26 shows a further decrease in average flux just before the transition 
into the B-mode, and perhaps truly turns off completely at these times. Furthermore, 
simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg 
telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the 
Q-mode and B-mode occurs within one single rotation of the neutron star, and that 
it is concurrent across the range of frequencies observed. 

Key words: stars: neutron - pulsars: magnetosphere - pulsars: individual: 

PSR B0823+26 - radio telescopes 


1 INTRODUCTION 

The average pulse profile shapes of radio pulsars are remark¬ 
ably stable when summed over at least a few hundred rota¬ 
tions. This reflects their clock-like rotational predictability 


(e.g. Liu et al.|2012) and enables their use as precision phys- 

ical probes (e.g. Weisberg & Taylor 2005 

Hobbs et al. 2010 

|Van Eck et al.|2011). Nonetheless, pulsar 

emission variabil- 


ity has been recognised on almost all timescales that are ob- 
servationally accessible (e.g. Keane] 2013): from nanosecond 


‘shots’, expected to be the quanta of pulsar emission (Han- 
kins et alT[2 003), to multidecadal variations presumably due 


to evolution of the magnetic field ( Lyne et al.|[ 2013). 

Indeed, there are a number of common emission phe¬ 
nomena which show that pulsar magnetospheres can be both 
stable and dynamic. For instance, drifting subpulses are a 
periodic drift of consecutive single pulses through the av¬ 
erage pulse window, at a rate of tenths to a few times the 
spin period ( Drake k, Craft|p~968 Weltevrede et al. 2006). 
Furthermore, hundreds of pulsars show nulling: a sudden 
cessation in emission between otherwise strong and steady 
pulses, spanning one to a few spin periods in most cases, but 


lasting upwards of several hours in some cases (e.g. Wang 
et al. 2007). Extreme nulling events can last for days to 


years (such sources are sometimes termed ‘intermittent pul¬ 
sars’ depending on the arbitrary length of time for which 
the pulsar is undetected). As with shorter nulls, the transi¬ 
tion between the detected emission, ‘ON’, and the null state, 
‘OFF’, occurs in less than 10 seconds, and perhaps within a 
single rotation (Kra mer et al.| 2006). For longer-term nulling 
pulsars, the spin-down rate during the ‘ON’ state has been 
detected to be approximately 1.5- to 2.5-times larger than 
that during the ‘OFF’ state (Kramer et al.||2006| [Camilo 


et al.||2012 Lorimer et al.||2012). 


Mode-changing is a similarly abrupt switch in the 
pulsed emission between two (or more) discrete and well- 


E-mail: sobey@astron.nl 


defined pulse profile morphologies ( Lyne|1971 Rankin et al. 
1988) that are observed to occur over the same time-scales 
as pulse nulling (i.e. seconds to years; e.g. Wang et al.| 2007 


Weltevrede et aQ|2011 Brook et al.||2014 ). This is also oc¬ 

casionally identified in variable linear ( Backer et al.|[l976 ), 
and circular polarisation ( Cordes et al.|1978 ). In some cases, 
pulse profile changes are also correlated with large changes 
in spin-down rates (0.3-13%) (Ly ne et aL||201 0). 

The above emission phenomena are thought to be con¬ 
nected, and related to changes in the current flows in the 


pulsar magnetosphere (Lyne et al. 2010 van Leeuwen & 


Timokhin 2012 Li et al. 2012b). For example, nulls can 


be regarded as a mode during which emission ceases or is 
not practically detectable. Most pulsars may exhibit these 
emission phenomena, but the effects of pulse-to-pulse vari¬ 
ability are hard to identify in sources with low flux den¬ 
sities, because summing over many pulses is required to 
achieve sufficient signal-to-noise ratio (S/N). Also, it is pos¬ 
sible that fewer pulsars are identified with longer character¬ 
istic timescale phenomena, such as mode-changing and ex¬ 
treme nulling, because frequent observations over a longer 
timespan are necessary to enable this. The hope is that bet¬ 
ter understanding these emission characteristics may lead 
to a better understanding of the magnetospheric emission 
mechanism itself. 

PSR B0823+26 is the focus of this work because it 
has been shown to exhibit a host of emission phenomena. 
This includes bursts of subpulses that drift slowly in longi¬ 
tude, towards both the leading and, more commonly, trailing 
edge of the main pulse (Weltevrede et al. p007| |. The nulling 
fraction has previously been estimated to be 7(2) and 
nulls occur clustered in groups more often than at random 
(Herfindal &; Rankin 2009 Redman fc Rankin||2009| . Re¬ 
cently, Young et al. 2012) showed that PSR B0823+26 nulls 


1 Note: Throughout the paper, numbers in parentheses are the 
uncertainties corresponding to the least significant figure in the 
value quoted. 
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over a broader range of timescales (from minutes to at least 
five hours), and that the pulsar is detected only 70-90% 
of the time when observed. They also found no detectable 
change in spin-down rate between emitting and long null 
states to a limit of less than 6% fractional change. Long¬ 
term timing of the pulsar shows relatively high timing noise 
residuals of 89ms ( Hobbs et al.||2 004), including changes 


the sign of the spin frequency second derivative (Shabanova 
et al.||2013 ). Thus, PSR B0823+26 is a good candidate for 


studying how the many facets of the pulsar emission mech¬ 
anism are observed and related in one source. However, ad¬ 
mittedly this may complicate matters for studying any one 
of these properties individually. 

PSR B0823+26 was one of the first pulsars to be dis¬ 
covered (Craft et al. 1968), and is one of the brightest radio 


400 MHz (Lorimer et al.|l995 

) and a distance of 0.32^|j^ kpc 

(Verbiest et al. 

2012 

). It has been detected from 40 MHz 

(Kuz’min et al. 

2007 

up to 32 GHz (Lohmer et al. 2008), 

with little profile evolution (Bartel et al.|1978). It is also one 


of ten slowly-spinning, rotation-powered pulsars with pulsa¬ 
tions detected in X-rays ( Becker et al.| 2004). The interpulse 
in the pulse profile suggests that it is nearly an orthogonal 
rotator, which is also supported by polarisation observations 
(Everett k Weisberg 2001). 

To further investigate the emission phenomena dis¬ 
played by PSR B0823+26, we conducted radio observations 
at low frequencies (<200 MHz) using the Low-Frequency 
Array (LOFAR), where brightness modulation due to scin¬ 
tillation is averaged out across the band. We also con¬ 
ducted simultaneous observations using the LOFAR, West- 
erbork, Lovell, and Effelsberg telescopes to record data from 
110 MHz to 2.7 GHz. 

The observations and data reduction are described in 
Section [2] In Section [3] we report on the data analysis and 
results of the observations. We show that the long-term nulls 
previously reported are in fact a very weak and sporadically 
emitting mode, similar to that identified in PSR B0826—34 
(Esamdin et al. 2005). Hereafter we refer to this mode as the 
‘quiet’ (Q) emission mode, and therefore the regular strong 
emission as ‘bright’ (B), analogous to other mode-changing 
pulsar^] In Section [I] we discuss these and further results 
and present our conclusions. 


2 OBSERVATIONS 

LOFAR observations of PSR B0823+26 were obtained on 
2011 November 13-14, 2012 February 9, and 2013 April 7 
(see van Haarlem et aL||2013 for a description of LOFAR). 
Between 5 and 21 High-Band-Antenna (HBA) core stations 
were coherently combined using the LOFAR Blue Gene/P 


2 PSR B0823+26 has a dispersion measure (DM) of 
19.454(4) pc cm -3 (Hobbs et al. 2004), and both diffractive 


and refractive scintillation is observed - the decorrelation band¬ 
width at 1.7 GHz is 81(3) MHz (Daszuta et al.[|| 2013| ). We note, 
however, that the 2-10 hour ‘disappearances’ of PSR B0823+26 
interpreted as due to scintillation by|D aszuta et af] ( |2013| ) are 
more likely to be instances of the weakly-emitting Q-mode 
observed in this work. 


correlator/beam-former to form a tied-array beam (see Stap- 
|pers et~ah ]|20TTJ for a description of LOFAR’s pulsar ob¬ 
serving modes). For a summary of the specifications of the 
observations, see Table [l] 

Starting at 23:18 UT on 2011 November 13, and ending 
at 08:12 UT on 2011 November 14, 26 three-minute obser¬ 
vations were taken with a gap of 19 minutes between suc¬ 
cessive observations. Data were taken using the central six 
Core stations (the ‘Superterp’), and were written as 32-bit 
complex values for the two orthogonal linear polarisations at 
a centre frequency of 143 MHz and bandwidth of 9.6 MHz. 
PSR B0823+26 switched to emitting in the B-mode during 
the thirteenth observation. The pulsar was not detected dur¬ 
ing the other 25 similar observations in this campaign, even 
after all observations were summed together in time. 

To further evaluate the moding behaviour at LOFAR 
frequencies, a longer, three-hour observation was taken on 

2012 February 9 using five of the central HBA Core stations. 
The linear polarisations were summed in quadrature and 
the signal intensities (Stokes I) were written out as 32-bit, 
245.76-ps samples, at a centre frequency of 143 MHz with 
47.6 MHz bandwidth. 

To investigate the broadband behaviour of 

PSR B0823+26, a further eight-hour observation on 

2013 April 7 was conducted simultaneously using LOFAR 
and three telescopes observing at higher frequencies: the 
Westerbork Synthesis Radio Telescope (WSRT), the Lovell 
Telescope and the Effelsberg 100-m Telescope, see Table [T| 
for a summary. The LOFAR observation used 21 HBA Core 
stations, where the linear polarisations were summed in 
quadrature, and the signal intensities (Stokes I) were writ¬ 
ten out as 32-bit 245.76-ps samples, at a centre frequency of 
149 MHz with 78 MHz bandwidth. The WSRT observation 
was taken in tied-array mode using the PuMall pulsar 


backend (Karuppusamy et al. 2008). Baseband data were 


recorded for 8 slightly overlapping 10-MHz bands, resulting 
in a total bandwidth of 71 MHz centred at 346 MHz. The 
Lovell Telescope observed the pulsar using the digital 
Reconhgurable Open Architecture Hardware (ROACH) 
backend. Dual polarisations were Nyquist sampled over a 
400-MHz band centred at 1532 MHz and digitised at 8 bits. 
The Effelsberg 100-m Telescope observation was conducted 
using the PSRIX backend (ROACH-board system for 
online coherent dedispersion) that recorded data at a centre 
frequency of 2635 MHz with 80 MHz bandwidth. 

The data from the LOFAR observations were con¬ 
verted to 8-bit samples offline. Data from all observations 
were coherently dedispersed, except the LOFAR observa¬ 
tions on 2012 February 9 and 2013 April 7 that were in¬ 
coherently dedispersed, and folded using the pulsar’s rota¬ 
tional ephemeris (Hobbs et al. 2004) and the dspsr pro¬ 
gram ( van Straten k Bailes| |2011), to produce single-pulse 
integrations. Radio frequency interference (RFI) in the pul¬ 
sar archive files was removed in affected frequency subbands 
and time subintegrations using both automated (paz) and 
interactive (pazi) programs from the P SRC HI VEj^] library 
(H otan et al.| 2004). The LOFAR polarisation data were cal¬ 
ibrated for parallactic angle and beam effects by applying a 
Jones matrix (see No utsos et al.|[20l5 for further details). 


3 see http://psrchive.sourceforge.net for more information. 
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The linear polarisation parameters in the LOFAR data from 
the 2011 November 14 observation were also corrected for 
Faraday rotation, see Section [3+] for further details. 


3 ANALYSIS AND RESULTS 

Here we describe the analysis and results of the observa¬ 
tions summarised in Table [l] PSR B0823+26 switched be¬ 
tween emission modes at least once during each observing 
campaign. Initially we focus on the time- and frequency- 
averaged properties at LOFAR frequencies, Section [3.1[ and 
at a range of frequencies, Section |3.2| In Section |T3| we 
focus on the single-pulse analysis of the data from the si¬ 
multaneous observations. 


3.1 Emission Modes Description 

Figure^ shows an overview of the emission characteristics of 
PSR B0823+26 over the course of the three-minute LOFAR 
observation taken on 2011 November 14, during which the 
pulsar switched to B-mode emission. This is the only ob¬ 
servation with recorded polarimetric information. Figure [l] 
panel d , clearly shows the sudden increase in total flux den¬ 
sity 1.69 minutes after the start of the observation. Panels e 
and / show that this is also the case for the polarised flux 
density. The B-mode continues for at least 160 single pulses 
until the end of the observation. This is an inadequate num¬ 
ber of pulses to obtain a stable average pulse profile, which is 
shown for comparison in Fig. [T] inset plots. The total dura¬ 
tion of the B-mode in this case is not known exactly. During 
all 12 preceding and 13 following three-minute observations 
no emission was detected, and we infer that PSR B0823+26 
was either nulling and/or emitting in the Q-mode. Assuming 
this was also the case during the 19-minute gaps between the 
observations, the duration of the Q-mode and/or nulling can 
be constrained to at least 4.4 h before and 4.75 h after the 
observation summarised in Fig. [T| This multi-hour duration 
is similar to the length of prolonged nulls also observed by 


Young et al. (2012). It is also possible that ‘flickers’ of emis¬ 


sion on the order of minutes in duration, such as that shown 
in Fig. [I] could have occurred in between observations. 

Inspection of the single pulses from Fig. ^ shows that 
the transition between the Q and B-modes is very rapid, oc¬ 
curring within one single rotational period. This is notable 
in the profiles of the single pulses and in the dramatic in¬ 
crease, by at least 60-times, in median S/N. We also note 
significant changes in the pulse width and morphology. Fig¬ 
ure ^ lower right , shows that the full width at half max¬ 
imum, iC 5 o, of the initial single pulses in the B-mode are 
very narrow (1.4°), compared to the average pulse profile at 
the same frequency (7°, see Fig. [l] inset plots , and below). 
Subsequent B-mode pulses increase in width and drift to¬ 
wards earlier and later pulse phases until an equally narrow 
pulse (1.1°) occurs 80 pulses later at 2.40min. This may be 
expected due to intrinsic variability, e.g., subpulse drifting. 
Moreover, the narrow pulses are located towards the trailing 
edge of the average pulse profile, similar to what is observed 
from PSR B1133+16 (Kra mer et al.|2 003). 

Further investigation of the polarisation properties was 
also carried out. The RM of the integrated pulse profile and 
the two narrow highly-polarised pulses was determined using 


RM-synthesis (Brentjens & de Bruyn 2005). The observed 
RM of the integrated B-mode pulse profile was determined 
to be 6.28(4) radm -2 . We find the RMs of the narrow single 
pulses to be consistent with this value. This RM value was 
used to correct the linear polarisation Stokes parameters for 
Faraday rotation in order to maximise the linear polarisa¬ 
tion in the average pulse profile. The RM due to the ISM 
alone was determined to be 5.25(8) rad m -2 , after correcting 
for ionospheric Faraday rotation (Sotomayor-Beltran et al.l 


2013). This is in good agreement with (within less than 2a), 


and more precise than, th e current ATNF p ulsar catalogue 
value of 5.9(3) rad m -2 (Manchester 1974). The DM and 
corrected RM values can be used to determine the Galactic 
magnetic field direction and magnitude parallel to the line- 
of-sight towards PSR B0823+26 (e.g. |Noutsos et ah] 2008), 
which was calculated to be 0.332(6) pG. These results are 
discussed further in Section 1+51 


Variable pulse-to-pulse fractional polarisation was also 
identified. The two narrow pulses previously identified show 
relatively high polarisation (especially circular) fractions, 
see Fig. [l] inset plots. The initial B-mode pulse is 55(5)% 
right-hand circularly polarised, and the narrow pulse at 
2.4min is 56(5)% left-hand circularly polarised. The initial 
B-mode pulse and subsequent narrow pulse are 10(5)% and 
38(5)% linearly polarised, respectively. The fractional circu¬ 
lar polarisation for these single-pulses is therefore consider¬ 
ably larger than for the average pulse profile, which is 6(5)% 
circularly polarised and 15(5)% linearly polarised. 

In order to further examine the moding behaviour, we 
observed PSR B0823+26 again for a continuous three-hour 
period using LOFAR on 2012 February 9. As shown in Fig. 
[2] lower panels , the pulsar switched from Q- to B-mode emis¬ 
sion after 2.71 h. The lower panels of Fig. [2] clearly demon¬ 
strate that weak emission is detected from the start to 2.5 h 
into the observation, at approximately the same location in 
pulse phase as that of the main pulse during the B-mode. 
This weak, but significant, Q-mode emission has not previ¬ 
ously been detected; this mode was previously thought to 
be a long-term null. Between the weak Q-mode at 2.55 h 
and the beginning of the B-mode at 2.71 h, the flux den¬ 
sity of the pulsar decreases even further and no emission is 
detected. PSR B0823+26 rapidly transitions to B-mode at 
2.71 h, again within one single rotational period, and con¬ 
tinues in this mode until the end of the observation. This 
is clear from the abrupt increase in S/N, see Fig. [5] lower 
right , and also in the pulse profile, see Fig. [2] upper panel. 
The ratio between the median S/N of single pulses in the B 
and Q-modes is approximately 170. 

Figure [2] upper panel , illustrates the average pulse pro¬ 
files of the Q and B-modes. The B-mode pulse profile shows 
the known components from the literature, i.e., a narrow 
main pulse (MP), postcursor (PC), and interpulse (IP, sep¬ 
arated by almost exactly half a rotation period from the 
MP) (see |Backer et al.|1973| ). The newly discovered Q-mode 
pulse profile is notably weaker. The location of the peak in 
flux density is also skewed towards slightly later pulse phases 
compared to the B-mode MP (+2(1)°, 2.9 ms). Table [ 2 ] in¬ 
cludes a summary of the relative intensity, width and loca¬ 
tions of the B and Q-mode pulse profile components. Details 
of how these were calculated are included in Section 13.21 
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Table 1 . Summary of observations. Columns 1-8 indicate the observing telescope, date, start time, integration time, centre frequency, 
bandwidth, individual channel width, and sampling time. Column 9 shows the duration of the B-mode relative to the start time, 
where ‘+’ indicates that the pulsar remained in the B-mode until at least the end of the observation. For the LOFAR observations on 
2011 November 14, 2012 February 9, and 2013 April 7: the observation IDs are L34789, L45754, L119505, respectively; HBA[0,1] stations 
included in tied-array mode were CS00[2-7], CS00[2,3,4,6,7] and CS0[01-07,11,17,21,24,26,28,30,32,101,103,201,301,302,401], respectively. 
Due to increased observing bandwidth and number of stations used, the 2012 February 9 and 2013 April 7 LOFAR observations were 2 
times and 10 times more sensitive than the original 2011 November 14 LOFAR observation, respectively. 


Telescope 

Date 

Start Time 

Tint 

Frequency Bandwidth Channel-width 

Tsamp 

Time in B mode 


(dd.mm.yyyy) (UT) 


(MHz) 

(MHz) 

(kHz) 

(H 


LOFAR 

14.11.2011 

03:42 

3.0 (min 

i) 143 

10 

12 

81.9 

1.69+ (min) 

LOFAR 

09.02.2012 

20:30 

3.0 (h) 

143 

48 

12 

245.8 

2.71+ (h) 

LOFAR 

07.04.2013 

14:00 

8.0 (h) 

149 

78 

12 

245.8 

3.07-3.11,5.09+ (h) 

WSRT 

07.04.2013 

14:00 

8.0 (h) 

346 

71 

24 

64.7 

3.07-3.11,5.09+ (h) 

Lovell 

07.04.2013 

14:00 

8.0 (h) 

1532 

400 

4000 

518.2 

3.07-3.11,5.09+ (h) 

Effelsberg 07.04.2013 

14:43 

7.3 (h) 

2635 

80 

195 

259.1 

(3.07-3.11),5.09+ (h) 



S/N _ Phase Phase _ Phase w 50 (deg) 


0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

Log[flux density] (normalised) 

Figure 1. Overview plot of the LOFAR observation on 2011 November 14. Panels a,b,c: Time- and frequency-averaged pulse profiles for 
the B-mode in total (Stokes /), linearly polarised (\L\ = Stokes y/Q 2 + U 2 ) and circularly polarised (Stokes V) intensities, respectively. 
The average profile before the start of the B-mode is not shown as it is indistinguishable from the noise in this short observation. Panels 
d,ej: The flux density of single pulses against pulse phase and integration time in total (increasing from white to black, see colour bars), 
linearly polarised (increasing from white to red) and absolute circularly polarised (increasing from white to blue) intensities, respectively. 
Lower left: S/N of total intensity single pulses against time. Lower right: Full-width-half-maximum, w$o, against time for total intensity 
single pulses with S/N>7. Inset plots: Normalised polarisation profiles for two narrow and highly polarised single pulses indicated 
with arrows: initial B-mode pulse (left), and second pulse with similar characteristics at 2.4min (right). Lower panels show the total 
(black), linearly polarised (red) and circularly polarised (blue) intensities. The stable average pulse profile in total intensity from the 
2012 February 9 LOFAR observation is also shown (dashed line) for the same frequency and using the same timing ephemeris (from Pilia 
et al. submitted). Upper panels show the position angle of the linear polarisation (P.A. = 0.5 tan ~ 1 (U/Q)). 
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LOFAR 143 MHz 



Figure 2. Overview plot of the LOFAR observation taken on 
2012 February 9. Lower left panel: Flux density of the pulsar 
against pulse phase and integration time (logarithmic scaling, in¬ 
creasing from white through blue to green, see colour bar at the 
upper right). The modes of emission are identified: Q-mode (0— 
2.55 h), B-mode (2.71+ h), and, in the period of approximately 10 
minutes in which the pulsar is not detected, tentatively ‘OFF?’ 
(2.55-2.71 h). These emission modes are delimited by the dashed 
lines. Subintegrations removed due to RFI are left blank, and the 
most affected time period is also labelled. Upper panel: Time- 
and frequency-averaged total intensity pulse profiles separated 
by emission mode into Q (0-2.55 h, grey line) and B (2.71+h, 
black line), normalised to the peak flux of the B-mode profile. 
The postcursor (PC) and interpulse (IP, 5x magnified) are also 
indicated. Lower right panel: S/N of 42-second total intensity 
sub-integrations against integration time. 


3.2 Multi-frequency Analysis 


In light of previous studies of other pulsars (e.g. Hermsen 
et al. 2013), mode-switching is expected to occur over a 


broad-range of wavelengths, even up to X-rays. To investi¬ 
gate whether this is the case for PSR B0823+26, we observed 
the pulsar using four telescopes (LOFAR, WSRT, Lovell, 
Effelsberg) at a range of centre frequencies (149, 345, 1532, 
2635 MHz) simultaneously. The overview of these eight-hour 
observations, Fig. [3J shows that both Q and B-modes were 
detected again, and that the behaviour of PSR B0823+26 is 
very similar across the range of frequencies. 

The lower panels in Fig.[3]show that from the beginning 
of the observations until approximately 2h, the weak Q- 
mode emission is visible at ah four frequencies. Between ~2 
and 3.07 h the Q-mode becomes even more weak and is prac¬ 
tically undetectable. After this, the pulsar abruptly ‘flickers’ 
into B mode emission for approximately 160 pulses. The pul¬ 
sar transitions back into the weakly-emitting Q-mode be¬ 
tween 3.11 and 4.8 h, which is most clearly detected at the 
two lower frequencies, but also towards the end of this pe¬ 


riod at the two higher frequencies. The flux in the Q-mode 
once again decreases even further and becomes practically 
undetectable between 4.8 and 5.09 h, just before the pulsar 
transitions again into B-mode for the remainder of the ob¬ 
servation. This is also demonstrated in terms of S/N against 
observing time in Fig. [4] 

Figure [3J upper panels , illustrate the B-mode profiles 
and the previously unknown and much weaker Q-mode pulse 
profiles between 110 MHz and 2.7 GHz. The B-mode MP, 
PC, and IP are visible at ah frequencies and, as expected 
from the literature, there is little profile evolution across 
the observed frequency range. Again, the Q-mode is much 
weaker compared to the MP of the B mode, and the average 
peak flux occurs towards slightly later pulse phases (> +1°, 
depending on frequency). 

On ah occasions where the emission switches from the Q 
to B-mode, the transition occurs within one single rotational 
period at ah frequencies, identified by the notable increase 
(70 x) in S/N of single pulses and change in pulse profile 
morphology. Figure [3] lower panels , also show that the mode 
switch from Q to B at 5.1 h occurs simultaneously for ah 
pulse profile components, i.e., the MP, PC and IP, at ah 
frequencies. 

To determine whether the mode-switch occurs coinci- 
dently across all frequencies observed, the MJD of the initial 
B-mode single-pulse after the mode change was determined 
for each frequency. The expected time delays of equivalent 
pulses due to dispersion in the cold plasma of the ISM be¬ 
tween the centre frequencies used at Effelsberg, Lovell, and 
WSRT and the centre frequency used at LOFAR are 3.63, 
3.61, and 2.97 s, respectively. Taking this into account, we 
find that the mode-switch is coincident at the frequencies 
observed to within 0.01s. 

On the single occasion where we observe the pulsar 
switch from the B to Q-mode, it is more difficult to pinpoint 
the transition because the flux density of single pulses de¬ 
creases gradually over a few rotational periods, similar to the 
behaviour during the 2011 November 14 observation shown 
in Fig. [I] panel d, rather than showing an abrupt change 
in S/N. For the three lower frequencies, B-mode emission 
was detected over a 160 single-pulse timescale during the 
B-mode ‘flicker’ at 3.1 h. Since the transition from the Q to 
B-mode was determined to be simultaneous across the ob¬ 
served band, it is probable that this is also the case for the 
B-to-Q-mode change. 

Figure [4] shows the S/N against time for the simulta¬ 
neous observations. This clearly demonstrates the simulta¬ 
neous wide-band transitions between emission modes, espe¬ 
cially the decrease in flux density in the Q-mode just before 
the rapid transition to the B-mode at 5.1 h. Here we again 
see that the B- and Q-mode differ in median flux density by 
over two orders-of-magnitude. 

Fig. [4] illustrates fluctuations in pulse total intensity 
due to scintillation, which are much more evident in the 
prolonged B-mode emission at the two higher observing 
frequencies. This is because the scintillation bandwidth 
becomes comparable to the observing bandwidth at the 
higher frequencies. The DM of PSR B0823+26, measured 
using the LOFAR observation on 2013 April 7 using the 
PSRCHIVE PDMP program (H otan et al.|2004 ), was found 
to be 19.475(3) pc cm -3 . This amount of dispersion con¬ 
tributes to PSR B0823+26 falling into the strong diffrac- 
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LOFAR 150 MHz 


WSRT 350 MHz 



0.2 0.3 0.4 0.5 0.6 0.7 

Log[flux density] (normalised) 
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Log[flux density] (normalised) 
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Figure 3. Overview of the multi-frequency simultaneous observations taken on 2013 April 7, using: LOFAR (upper left), WSRT (upper 
right), Lovell (lower left), and Effelsberg 100-m (lower right). Upper panels: Time- and frequency-averaged total intensity pulse profiles 
separated into Q-mode (0-3.07h, 3.11-5.09h; 10X magnified; grey lines) and B-mode (3.07-3.11 h, 5.09+h, black lines), normalised to 
the peak flux of the B-mode profile. There also appear to be periods in which the pulsar is not detected (‘OFF?’; ~2-3.07h, 4.8-5.09 h). 
The PC (lOx magnified in the lower plots) and IP (10X magnified) are also indicated. Lower panels: Flux density of the pulsar against 
pulse phase and integration time (increasing from white through blue to green, see colour bars). Subintegrations removed due to RFI 
are shown as white. The modes of emission identified (Q-, B-, and, tentatively, ‘OFF?’) are also labelled, and these are delimited by 
the dashed lines. The off-pulse noise level in the LOFAR observation shows an increase around the transition from Q- to B- to Q-mode 
emission around 3h because there was a period of increased RFI around this time, and more conservative RFI flagging was done to 
conserve as much data around this period as possible. Modulation of the B-mode brightness due to scintillation is also visible in the 

higher frequency Lovell and Effelsberg data. 
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Figure 4. S/N of total intensity pulses versus integration 
time from the multi-frequency simultaneous observations of 
PSR B0823+26. Upper to lower panels: Effelsberg 100-m at 
2635 MHz using 80-s subintegrations; Lovell at 1532 MHz using 
140-s subintegrations; WSRT at 345 MHz using 140-s subinte¬ 
grations; LOFAR at 149 MHz using 110-s integrations. Discon¬ 
tinuation of lines indicates flagged subintegrations due to RFI. 
The periods of B-mode emission are indicated by the shaded grey 
backgrounds and labelled ‘B\ 


tive scintillation regime when observed below approximately 
5 GHz (Daszuta et al. 2013). At lower frequencies the scin¬ 
tillation bandwidth is much smaller, causing its effects to be 
averaged out more effectively across the band. 

The B- and Q-mode multi-frequency average pulse pro¬ 
files obtained from the 2012 November 9 and 2013 April 7 
observations, see Figs.[2]an di were used to compare the be¬ 
haviour of the pulse profile components with frequency. The 
profile components were fitted with Gaussian functions us¬ 
ing the Levenberg-Marquardt algorithm (Levenberg 1944) to 
provide the normalised peak amplitude, peak location, icso, 
and full width at tenth of maximum, w io, values, and their 
respective errors. Table [2] provides the summary of these 
characteristics for the pulse profile components at multiple 
frequencies for B and Q-mode emission. Peak amplitude is 
measured in percentage relative to the B-mode MP. 

After quantifying characteristics of the profile features 
at each frequency listed in Table [2] the trend with respect 
to frequency was explored. Assuming that each feature, x , 
follows a power law with respect to frequency, x oc z/ 77 , a 
weighted best fit was performed to determine the appropri¬ 
ate value for the exponent, r\. The resulting values and errors 
for the exponents are shown in the bottom row of Table [2] 

We find that the width of all pulse profile components 
decrease with increasing frequency in both emission modes, 
although this trend is most notable for the Q-mode pulse. 
This is consistent with previous results for the B-mode pulse 
profile (e.g. |Kramer|1994 ), and also with radius-to-frequency 
mapping (RFM; e.g. Gil et al. 2002). The largest intra¬ 
channel DM smearing at the lowest LOFAR frequency ac¬ 
counts for 0.3° in pulse longitude, which does not affect the 
conclusions obtained from the low-frequency pulse profile. 
The B-mode PC and IP move further from the MP and 



Figure 5. Low-frequency spectra of PSR B0823+26, separated 
into B- (filled circles) and Q-mode emission (open squares), ob¬ 
tained from the LOFAR observation taken on 2013 April 7. 
Frequency-axis error bars show the bandwidth used to calculate 
each point. The dashed lines show the best power-law fit to the 
data. 


decrease in amplitude relative to the MP with increasing 
frequency. The positions of the PC and IP and their fre¬ 
quency dependence are in agreement with previous results 
(Hankins & Fowler 1986). This is further discussed below in 
Section [4] Table [2] shows that the amplitude of the PC has 
the strongest frequency dependence and the location of the 
IP is the least dependent on frequency. 

The relative pulse amplitude of the Q-mode profile is 
approximately equal to the amplitude of the B-mode IP, 
and also decreases with increasing frequency. The peak am¬ 
plitude of the Q-mode is located towards later pulse phases 
compared to the B-mode MP at the same frequencies, and 
seems to arrive progressively earlier in pulse phase towards 
higher frequencies. Furthermore, Table [2] shows that both 
B- and Q-mode pulse profiles show no significant change 
over one year between the LOFAR observations from 2012 
February 9 and 2013 April 7. 

To further investigate the relative properties of the B 
and Q-modes at low-frequency, the spectrum in each mode 
was calculated using the eight-hour LOFAR observation. 
The flux density, S, as a function of frequency was obtained 
using the radiometer equation, including an empirically- 
derived and tested formula for the effective area of the LO¬ 
FAR stations, and correction for the primary beam using 
the theoretically-derived and tested beam model, described 
in Hassall, et al. (in prep.), and references therein. The B- 
mode and Q-mode data were divided into 5 and 3 subbands, 
respectively, deemed to be appropriate to optimise the re¬ 
spective S/N. A line of best fit was calculated assuming a 
single power law across the bandwidth of the observation, 
S oc z/*, and the error was determined considering the flux 
density errors only. Figure [5] shows the resulting spectra in 
the two emission modes. 

The difference in flux density between the B and Q- 
modes is striking, see Fig. [5] In fact, the average flux density 
in the B-mode, 370 mJy at 149 MHz, is almost two orders-of- 
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Table 2. Summary of the pulse profile components from the LOFAR observation on 2012 November 9 (LOFARi) and the simultaneous 
multi-frequency observations on 2013 April 7. Columns 1 and 2 indicate the telescope and central frequency used. 1^50 of the B-mode MP, 
PC, and IP, and the Q-mode pulse are shown in columns 3, 7, 11, and 14, respectively. w\q of the B-mode MP and PC are summarised 
in columns 4, and 8 , respectively. The peak flux location of the B-mode PC and IP, and the Q-mode pulse, with respect to that of the 
B-mode MP (where £ +’ indicates later pulse phases) are shown in columns 5, 9, and 12, respectively. The relative amplitude of the peak 
flux of the B-mode PC and IP, and Q-mode pulse, with respect to the MP (corrected for different integration times using Table Q are 
summarised in columns 6 , 10, and 13, respectively. The bottom row shows the frequency dependence of each of the profile features as 
the exponent of a power law, 77 . 




B-mode MP 


B-mode PC 

B-mode IP 



Q-mode 


Telescope 

Frequency 1050 

m 0 

location 

amplitude W 50 

w 10 location 

amplitude w$o 

location amplitude 

W 50 

(Name) 

(MHz) 

(°) 

(°) 

(°) 

(%) 

(°) 

(°) (°) 

(%) 

(°) 

(°) 

(%) 

(°) 

LOFARi 

143 

6.68(5) 

14.06(8) +31.5(1) 

3.76(5) 

9.8(3) 

17.9(5) -178.5(7) 

0.90(5) 

20.5(6) 

+ 2 ( 1 ) 

0.89(5) 

9.4(2) 

LOFAR 

149 

6.73(3) 

13.36(8) +31.5(1) 

3.19(5) 

9.6(2) 

17.4(2) -179.3(2) 

0.75(5) 

16.2(5) 

+1.9(2) 

0.62(5) 

9.6(4) 

WSRT 

346 

4.29(3) 

9.14(8) 

+31.7(1) 

2.09(5) 

9.4(2) 

17.2(3) -180.0(3) 

0.30(5) 

15.2(2) 

+ 1 . 8 ( 2 ) 

0.47(3) 

8.7(5) 

Lovell 

1532 

3.12(3) 

6.33(8) 

+36.0(2) 

0.69(4) 

8.3(2) 

15.2(2) -181(1) 

0.23(6) 

12.3(8) 

+1.3(5) 

0.39(9) 

3.7(3) 

Effelsberg 2635 

2.89(7) 

5.97(9) 

+36.1(1) 

0.51(3) 

7.2(3) 

13.2(3) -182(1) 

0.12(3) 

11 ( 1 ) 

+0.9(6) 

0.7(3) 

1.8(4) 

All 

77 

-0.33(4) 

-0.30(4) +0.050(9) 

- 0 . 66 ( 6 ) 

-0.08(2) 

-0.08(2) +0.005(1) 

-0.63(2) 

-0.19(6) 

-0.16(4) 

-0.3(2) 

-0.37(9) 


magnitude greater than that of the Q-mode, 5 mJy. This is 
probably a factor in the previous non-detection of this weak 
Q-mode. 

The exponent, <a, obtained from the continuous power- 
law best fits, weighted by the uncertainties, for the B and 
Q-modes are -2.6(5) and -2.3(1), respectively. Therefore, 
there is no significant difference in the power-law index of 
the spectra between emission modes. This is further dis¬ 
cussed in Section [4] The relatively steep spectral index and 
strong scintillation are also possible factors in the previous 
non-detection of the weak Q-mode emission at higher fre¬ 
quencies. The increased bandwidth used for the observa¬ 
tion using the Lovell Telescope on 2013 April 7 in this work 
(400 MHz) also provided almost a factor of two better sen¬ 
sitivity in comparison to the observations in | Young et al.| 

p0l2| . 


3.3 Single-pulse analysis 

Potentially interesting emission features were investigated 
and compared for the B and Q-modes by analysing sin¬ 
gle pulses from the eight-hour simultaneous observations. 
This was conducted using a software package provided by 
P. Weltevrede (private communication) to construct cen¬ 
tred, gated, de-baselined, and further RFI-excised single¬ 
pulse archives (e.g. | Weltevrede et ah||2006| ). This enabled 
us to extract information regarding nulling, subpulse drift¬ 
ing, and pulse energy distributions (PEDs). 

3.3.1 Pulse Energy Distributions 

Longitude-resolved PEDs were calculated for each pulse pro¬ 
file component of the B and Q-modes by integrating over a 
selected pulse longitude bin range, see Fig. [6] Pulse longitude 
ranges corresponding to w\q were selected to optimise S/N. 
These are indicated underneath the average pulse profiles 
in Fig. [6] upper panels. The longitude range of the off-pulse 
(OP) region shown in Fig. [6] was selected to correspond to 
the range selected for the Q-mode pulse, and to be located 
mid-way between the B-mode MP and IP. 


Figure [6] upper panels , shows the normalised flux- 
density profiles of the B and Q-mode. Note the difference 
in the off-pulse noise level. The lower panels in Fig. [6] il¬ 
lustrate the normalised PED with respect to the average 
pulse energy, <E>, shown separately for the Q (left) and 
B-modes (right). The B-mode MP and PC PEDs are clearly 
distinct from that of the off-pulse region and resemble log¬ 
normal probability density distributions, see inset panel. Us¬ 
ing a method similar to that in Weltevrede et al. (2006), the 
best lognormal fit to the PED for the MP for pulse ener¬ 
gies 1.5-20 times greater than the average yields a standard 
deviation of 0.77(15). The IP PED is comparable to that 
found for the off-pulse region because of its relatively weak 
single-pulse flux. 

Q-mode emission shows a distinctly different PED, 
which falls off much more rapidly at low energies and re¬ 
sembles a power-law distribution. The Q-mode PED is less 
distinguishable from that of the off-pulse region at lower 
pulse energies due to the much weaker emission, but shows 
a tail of less frequent higher-energy pulses. The best power- 
law fit to the PED for pulse energies 1-8-times greater than 
the average gives an exponent of-3.2(2). 

To investigate how the energy of the pulses affects the 
average pulse profile, we constructed integrated PEDs from 
the longitude-resolved PEDs, shown in Fig. [6] to obtain the 
average pulse profile for a number of pulse energy ranges. 
We find that the most energetic B-mode single pulse profiles 
show PCs located towards lagging pulse phases. Moreover, 
these pulses are approximately eight-times brighter com¬ 
pared to the average pulse profile. Despite this, the aver¬ 
age pulse profile shows remarkably little pulse profile evolu¬ 
tion with increasing pulse energy, similar to that observed 
in PSR B1133+16 (Kramer et al. 2003). In addition, sin¬ 
gle pulses with the highest flux in the Q-mode have smaller 
W 50 values than the average pulse profile. There is also a 
more notable contrast in pulse profile evolution with energy 
in the Q-mode. That is, pulses with energies greater than 
twice the average possess profiles akin to a skewed Gaus¬ 
sian, similar to that in Fig. [6] top right. Pulses with energies 
less than twice the average have profiles that are weak (dis- 
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Figure 6. Upper panels: Normalised average pulse profiles (black line) for the Q-mode (left) and B-mode (right), showing longitude 
ranges included in the pulse-energy analysis (delimited by the black dotted lines). Lower panels: Normalised longitude-resolved PEDs 
for the Q- (left) and B-mode emission (right), using the LOFAR observation taken on 2013 April 7. The MP, PC, IP, and off-pulse (OP) 
distributions are shown (red, blue, green, and black dotted lines, respectively). Inset panel: B-mode MP PED data (red points) and best 
lognormal fit to the data (red line), plotted on logarithmic axes. 


cernible from the baseline at approximately 3-cr), irregular, 
and somewhat top-hat-like in shape, for examples see Fig. 
[7j left. Hence, the pulse profiles in the Q-mode seem to show 
much greater variability with energy compared to those in 
the B-mode. 


3.3.2 Nulling 

The integrated PEDs were used for complimentary analy¬ 
sis of the nulling fraction (NF). We estimated the NFs by 
subtracting the off-pulse pulse-energy samples from the on- 


pulse PEDs, as described in Wang et al. (2007). The analysis 


was undertaken separately for the Q-mode, the 160-pulse B- 
mode ‘dicker’, and the prolonged B-mode from the 7 April 
2013 LOFAR observation, summarised in Fig. [3] 

PSR B0823+26 displays nulling in both emission 
modes. More specifically, the NF for the prolonged B-mode 
was determined to be 1.8(5)%. This is lower than, but still 
within 2.2-cr of, previous estimates (Herhndal & Rankin 
|2009| . The NF for the B-mode ‘flicker’ was found to be 
15(1)% - almost ten times larger than that during the pro¬ 
longed B-mode emission. 

Further inspection of the single pulses during this short¬ 
lived period of B-mode emission also indicates that the NF 
increases towards the transition back to the Q-mode. Inspec¬ 
tion of simultaneous pulses from the 346 MHz and 1.5 GHz 
data also indicate that the nulls occur simultaneously at all 
frequencies, which is consistent with mode-switching, and is 
observed in several other cases (e.g. Biggs 1992; Hermsen 
et al.|2013 ). The NF throughout both occurrences of the Q- 


mode during the observation was determined to be 80(9)%. 
This indicates that there are over forty-times more nulls 
during the Q-mode than during B-mode emission. However, 
we note that the apparent nulls (in both B- and Q-mode 
emission) may be yet another weaker emission state that is 
currently undetectable. 

Figure [7] left, demonstrates the nature of the Q-mode 


emission through a 250-single-pulse extract from 4.6 h into 
the observation. A pulse extract from the B-mode emission 
is also shown in Fig. [7] right, for comparative purposes. 
The weak Q-mode emission consists of sporadic, mainly 
low-intensity bursts which occur in groups of approximately 
one to five pulses with no immediately identifiable periodic¬ 
ity. Therefore PSR B0823+26 seems to display similar be¬ 
haviour to the weak-mode emission identified in PSR B0826- 


34 (Esamdin et al. 2005), also somewhat resembling emission 
from rotating radio transients (RRATs) (McLaughlin et al. 
2006). This is also further discussed in Section |4| 


3.3.3 Subpulse Modulation Features 

The subpulse modulation properties of the prolonged B- 
mode at low frequency were investigated using the LOFAR 
observation on 2013 April 7, and the method described in 


Weltevrede et al. (2006). The aim was to compare the mod¬ 


ulation patterns exhibited by the different modes (e.g. _ 

Leeuwen & Timokhin 2012). Pulse stacks, for example see 


Fig. [7| were used to obtain the Longitude Resolved Fluctu¬ 
ation Spectrum (LRFS; [Back er 1970) and hence the vertical 
separation between drift bands, P3, measured in pulse pe¬ 
riods, Po. The LRFS obtained from the B-mode emission, 
Fig. [8] right, shows two drift features that are broader than 
0.05 cycles per period (CPP). One feature is located at the 
alias border P0/P3 = 0, and the other feature relates to P3 
— 5.26(5)Po- Since P3 is generally independent of observing 
frequency, this is in excellent agreement with the value of 


5.3(l)Po previously determined at 92-cm (Weltevrede et al. 

[2M7). 

We also constructed Two-Dimensional Fluctuation 
Spectra (2DFS; [Edwards & Stappers 2002) to determine 
the horizontal separation between the drift bands, P 2 , mea¬ 
sured in pulse longitude. For the B-mode emission we find 
that although the feature in the 2DFS peaks at zero CPP, 
the wings are asymmetric, indicating that drifting subpulses 
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Q-mode 150 MHz 


B-mode 150 MHz 



0.65 0.70 0.75 

Phase 



0.65 0.70 0.75 0.80 

Phase 


Figure 7. Representative examples of single-pulse profiles from 
the LOFAR observation on 2013 April 7, arbitrarily normalised. 
Left: sporadic and weak emission during the Q-mode, right: 13- 
mode regular emission for comparison. Single pulses rejected due 
to RFI are not shown. 


travel towards later pulse phases more often. Calculating the 
centroid of the asymmetric wings at 20% of the maximum 
allowed us to extract the value of P 2 — +90(20)°. This is in 
good agreement with the value obtained at 92-cm, +70° 


Weltevrede et al. (2007). Unlike P 3 , P 2 is found to vary with 
observing frequenc y ([E dwards &^tap£ers 2003). Since P 2 = 


+55°+^° at 21-cm (Weltevrede et al. 2006), the value deter¬ 


mined here also seems to be in agreement with the general 
trend that P 2 decreases with increasing observing frequency 
for PSR B0823+26. 

The same subpulse modulation analysis was also con¬ 
ducted for the prolonged Q-mode emission from the 2013 
April 7 LOFAR observation. No significant feature was de¬ 
tected in either the LRFS, Fig. [8] left , or 2DFS, indicat¬ 
ing that there is variability on all fluctuation frequencies. 
This may be attributed to the greater nulling fraction and 
very weak single pulses. However, PSR B0943+10, which 
shares many emission characteristics with PSR B0823+26, 
also shows subpulse drifting during the B-mode and not in 


Q-mode, which is also comparatively ‘disordered’ (Backus 
et al.|201 1). 


The longitude-resolved modulation index, a measure of 
the factor by which the intensity varies from pulse to pulse, 
was derived from the LRFS (Weltevrede et al. 2006), see 


Fig [8] upper panels. The median modulation index of the B- 
mode MP was found to be 1.4(5). The median modulation 
index of the PC component was found to be slightly larger, 
1.6(9). The B-mode MP and PC components show similar 
longitude-resolved modulation characteristics, although the 
PC shows slightly more variability than the MP. The me¬ 
dian modulation index of the Q-mode pulse was found to be 


over ten-times larger than the B-mode MP, 17(5), and shows 
somewhat different longitude-resolved modulation compared 
to the B-mode, see Figi upper panels. 

The pulse profile stability during B-mode was inves¬ 
tigated using the 7 April 2013 simultaneous LOFAR and 
WSRT data. The Pearson product-moment correlation co¬ 
efficient, p, was calculated between an analytic profile of 
PSR B0823+26 and the pulse profile integrated from an in¬ 
creasing number of effective single pulse numbers, A^ e f c , as 


described in Liu et al. (2012). We find that towards higher 


effective pulse numbers (>200-500) the LOFAR and WSRT 
data follow the expected trend, i.e., 1 —p oc N~f. (Liu et al. 
2011). However, there are some notable deviations from this 
trend, which tend to coincide with narrow, bright pulses. 
This suggests that over many pulse periods the pulse profile 
is stable, but over short periods the emission shows intrinsic 
pulse-to-pulse shape variations. This is not unexpected due 
to previous results that show the array of emission phenom¬ 
ena exhibited by PSR B0823+26. 


4 DISCUSSION AND CONCLUSIONS 

In this work we have confirmed that PSR B0823+26 shows 
a host of emission characteristics over a wide range of 
timescales. The most surprising finding from the work pre¬ 
sented here is that the long-term nulls, recently found by 


Young et al. (2012), are in fact a very weak and sporadi¬ 


cally emitting mode, which we refer to as Q-mode. We find 
that this mode is over 100 times weaker than the B-mode. 
We also find evidence for a further decrease in flux just be¬ 
fore the switch to B-mode, see Figs. [2] and [3] It is possi¬ 
ble that this could be a third emission mode, whereby the 
flux density of emission decreases even further or emission 
completely ceases. It would also be quite interesting if this 
weakest (possibly completely off) mode is found to always 
occur before the start of B-mode - as it has in the three 
instances we have observed. 

Considering the many studies which have been con¬ 
ducted on PSR B0823+26, it is surprising that the long-term 
nulling was only recently discovered ( Young et al.|| 2012). It 
may be that the emission from PSR B0823+26 has changed 
over a multi-year timescale, such that mode changing is more 
frequent or longer in duration than in previous years, though 
that remains to be shown. If so, this may be due to a multi- 
decadal change in emission, perhaps related to that observed 


in the Crab pulsar due to magnetic field evolution (Lyne 
e t al.|2013) . 


4.1 Magnetospheric Switches 

Upon investigation of the newly-discovered Q-mode pulse 
profile at frequencies between 110 MHz and 2.7 GHz, we 
find the pulse profile comprises a single component located 
within the regular B-mode MP envelope, but with a peak lo¬ 
cated towards slightly later pulse longitude, and a lower flux 
density by approximately two orders of magnitude. The PC 
and IP that are present during the B-mode are not detected 
during the Q-mode. A converse example of a pulsar that ex¬ 
hibits mode-changing, PSR B0943+10, has a pre-cursor in 
the Q-mode and not in the B-mode at 320 MHz (Hermsen 
et al.|2013 ). 
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Figure 8 . Pulse variability analysis products from the LOFAR observation on 2013 April 7, for the Q-mode (left) and the B-mode 
(right). Upper panels: the normalised integrated pulse profile (solid line), the longitude-resolved standard deviation (open circles), and 
the longitude-resolved modulation index (dimensionless quantity shown as the solid line with error bars). Lower panels: the LRFS is 
shown, with its horizontal integration to the left. 


During the Q-mode, the emission is weak and sporadic, 
with a very high NF and modulation index. This behaviour 
seems similar to the Q-mode identified in PSR B0834-26 
(Esamdi n et al.| |2005). It also seems reminiscent of RRAT- 
like emission that may imply a common, or at least re¬ 
lated, emission mechanism. This has previously been pro¬ 
posed from both theoretical (e.g. Jones 2013) and observa¬ 


tional perspectives (e.g. Burke-Spolaor & Bailes 2010; 


Keane 


et al.poTT ). The detection of the weakly-emitting Q-mode 


from PSR B0823+26 could mean that other pulsars cur¬ 
rently identified as extreme-nullers may also show similar 
weak emission if higher-sensitivity observations can be done 


(Esa mdin et al.|2005 Wang et aL||2007 Young et al.|2014 ). 

During the simultaneous observation on 2013 April 
7, the initial ‘OFF?’ period lasted approximately 0.87h, 
and preceded the short-duration 160 single-pulse B-mode 
‘flicker’. Later in the same observation, and also in the ob¬ 
servation on 2012 February 9, the ‘OFF?’ period was much 
shorter in duration, and the subsequent B-mode emission 
was much longer in duration. This further decrease in radio 
emission, which tentatively may be considered as a further 
separate emission mode, could be caused by a change in the 
particle density within the magnetosphere that soon after 
causes the switch to B-mode. We also speculate that the 
duration of the undetected emission towards the end of the 
Q-mode, just before the mode-switch to the B-mode, may be 
inversely correlated with the length of time the pulsar emits 
in the subsequent B-mode. We stress that this is speculative 
due to the limited sample, and more work will be done in 
future to investigate this by using longer and more frequent 
low-frequency observations. 

Single-pulse analysis allowed us to determine that the 
transition from the Q-mode to the B-mode occurs within 
one rotational period, and that it is broadband over the 


radio frequencies observed, similar to other mode changing 
pulsars (e.g. PSR B0943+10, Hermsen et al. 2013| Bilous 
et al.||2014|, a nd nulling pulsars (e.g. PSR B0809+74, Bar- 
tel et al.|1981 ). Although, considering the tentative ‘OFF?’ 
mode just before the switch to B-mode emission mentioned 
previously, the underlying mode change mechanism may be 
more complex than is evident from the single-pulse data. 
Within the theoretical framework of force-free magneto¬ 
spheres, it has been shown that multiple stable solutions 
with different structures are possible, e.g., with different 


sizes of closed field line region (Contopoulos 2005 Timo- 


khin||2010 ). Therefore, the B and Q-modes could be caused 


by two stable magnetospheric states with different magneto- 
spheric structures, affecting the current of relativistic parti¬ 
cles, and therefore the broadband radio emission mechanism. 
However, the mechanism that generates the wide range of 
observed timescales (e.g. each emission state lasting many 
hours and mode-switches occurring within one rotational pe¬ 
riod) is still elusive. 

4.2 Considering Spin-down 

PSR B0823+26 shows no detectable change in spin-down 


rate between modes (Young et al. 12012), similar to other 


pulsars which exhibit mode-changing on similar timescales 
e.g., PSRs B0834-26 and B0943+10, whilst other extreme¬ 
nulling pulsars identified as ‘intermittent’ often do (Kramer 
et al.|2006 Camilo et al.|2012 Lorimer et al.|20l2 ). In part, 
this is because it is difficult to detect small changes in spin- 
down if transitions occur more frequently than for the more 
extreme nulling pulsars (Young et al.| 2012|. This may also 
be due to the electromagnetic spin-down torque exert ed on 
the neutron star, proportional to sin 2 0 + (1 — k)cos 2 0 (Con- 
topoulos &; Spitkovsky|[2006 ), where 6 is the alignment an- 
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gle between the spin axis and magnetic moment, and k 
the ratio between the angular velocity of the pulsar and 
the angular velocity at the death line. Namely, if there is 
some small change in 6 between emission modes, then the 
change in spin-down torque is less severe in cases of nearly 


orthogonal rotators (e.g. PSR B0823+26, Everett V Weis- 


berg 2001). This is less certain in cases of nearly aligned 


rotators (e.g. PSR B0943+10, Deshpande & Rankin 2001) 
(Li et al. 2012a). 


4.3 Emission Characteristics 


Although we find that the spectral indices of the B and 
Q-modes are consistent, within errors, the broad-band spec¬ 
trum of PSR B0823+26 shown in Hassall et al. (in prep.) in¬ 
dicates a possible spectral turn-over at 127(25) MHz. Future 
observations using the LOFAR Low Band Antennas (30- 
90 MHz) would be ideal for further investigating the spec¬ 
tral properties of PSR B0823+26 in both emission modes (if 
the Q-mode is detectable), including providing constraints 
on the presence of a spectral turn-over, and hence supply¬ 
ing more clues about the origin of the difference between 
emission modes. 

PSR B0823+26 also shows certain emission characteris¬ 
tics which are reminiscent of several other pulsars. The PCs 
associated with the highest energy single pulses that are 
8-times brighter than average are reminiscent of the Vela 
pulsar’s (PSR B0833—45) ‘bump’ component after the main 
pulse, where ‘giant micro-pulses’ with flux densities exceed¬ 


ing 10-times that of the mean were observed (Kramer et al. 
2002). For comparison to the highly polarised single pulses 


detected in the 2011 November 14 LOFAR observation, PSR 
B0656+14, which may be a nearby RRAT, also shows strong 
highly-polarised single pulses limited to the leading and cen¬ 


tral regions of the average pulse profile (Weltevrede et al. 


2006). Moreover, narrow pulses observed after the transition 


to the B-mode occur towards the trailing edge of the aver¬ 
age pulse profile, similar to that observed in PSR B1133+16 



Polarisation observations have allowed the MP to be 
classified as a core component and the IP to be classified 
as a possible core component of the opposite magnetic pole 


(Weisberg et al. 1999 and references therein). The PC’s high 
linear polarisation fraction could originate from induced 
scattering of the main pulse emission into the background 
by the particles of the ultra-relativistic highly magnetised 
plasma (Petrova 2008). This could also explain the trend in 
the PC location growing more distant with increasing fre¬ 
quency, which is contradictory to the behaviour expected 
from RFM (e.g. Mitra V Rankin||2002 ) and has lead to re¬ 
luctance in its classification as a conal component (Weisberg 
et al.|[l999 ). We note that the decrease in the pulse profile 
component widths with increasing frequency is consistent 
with the RFM model (e.g. Gil et al. 2002^, although this 
may also be caused by other mechanisms (e.g. birefringence, 
McKinno n|[T9 97). 

It seems that the emission behaviour observed from 


PSR B0823+26 is very diverse, and may prove one of the 
most difficult pulsars to model in terms of accounting for 
the host of emission characteristics that arise over many dif¬ 
ferent timescales. However, observations of PSR B0823+26 
may also help in efforts to unify the diverse phenomenology 
of pulsar emission and provide vital data in terms of studies 
of the elusive pulsar emission mechanism. 


4.4 Pulsar Population Implication 

In the period-period-derivative (P-P) diagram for non- 
recycled pulsars from the ATNF pulsar catalogue (Manch- 
ester et al.|2005 ), PSR B0823+26 is located near the centre 
of the distribution, and therefore appears to be a very regu¬ 
lar pulsar in this parameter space. Other pulsars which have 
been identified to display at least one of the numerous emis¬ 
sion characteristics investigated here are located throughout 
most of the parameter space, apart from the area occupied 
by very young pulsars. Therefore, there is no obvious corre¬ 
lation with surface magnetic field strength or characteristic 
age (e.g. Wel tevrede et al.][2 007). Most pulsars may exhibit 
these emission phenomena, but the effects of pulse-to-pulse 
variability are diminished due to summing over many pulses, 
which is required to achieve sufficient S/N in the case of 
sources with low flux densities. 

To date, only a handful out of the ~ 2300 known pul¬ 
sars are identified as intermittent or extreme millers. This 
is because this requires frequent observations over several 
months/years to identify this long-term behaviour. This is 
also because they may also be difficult to identify in a pulsar 
survey due to the nature of the weak/undetectable emis¬ 
sion during quiet modes. There may also be pulsars that 
emit in a Q-mode-like state over longer timescales or con¬ 
tinuously, for which more sensitive radio-telescopes will be 


needed, such as the Square Kilometre Array (SKA, Garrett 
et al. 2010). Studying pulsars that display many of these 


emission characteristics, such as PSR B0823+26, will pro¬ 
vide further insights into the possible radio emission charac¬ 
teristics, the relationship between phenomena, and under¬ 
standing of the magnetospheric emission mechanism from 
pulsars and RRATs. 


4.5 Using Pulsars as ISM Probes 


Using the LOFAR observation on 2011 November 14, we 
found that there is no significant variation in the observed 
RM over small timescales due to changes in the pulsar mag¬ 
netosphere. This is reassuring in terms of using pulsar RMs 


as probes of the Galactic magnetic field (e.g. Van Eck et al. 
2011). Although the pulsar magnetosphere is assumed to 


have a minor contribution to the observed RM (Noutsos 
et~aT| 2009), an interesting area of future investigation may 
be whether the RM, or polarisation angle, as a function of 
pulse phase is in any way dependent on the emission mode. 
LOFAR’s low frequency and large fractional bandwidth can 
facilitate such studies, where high RM precision is impor¬ 
tant. Studies of the emission characteristics from pulsars, 
including in polarisation, will also continue to aid in better 
understanding the pulsar emission mechanism itself, so that 
pulsars can be used increasingly effectively as probes of the 
ISM. 
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